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Abstract 

The observation of several low energy events during SN1987A burst has an important 
weight in the attempts to learn about the properties of the supernova neutrinos, but is 
somewhat puzzling if compared with the expectations. In this work, we study the low 
energy events observed by Kamiokande-II, and consider in particular the possibility that a 
few of them are due to background. We focus our attention not only on the event 6, that 
falls below the energy threshold of 7.5 MeV, but also on the other four events that at 1 a fall 
below the energy threshold, namely the events 3,4,10,12. The a priori expectations on the 
number of background events above threshold would suggest 1 or 2 background events only. 
The volume distribution of the Kamiokande-II events is not a uniform distribution at 3 a, 
that suggests the presence of background events close or at the border of the volume used for 
the analysis, including the events 3,4,10. Next we checked the expected energy distribution 
assuming that the signal is due to v e p — > e + n and that the average antineutrino energy is 14 
MeV. The agreement with the observations is not perfect, but it is acceptable at the 11 % 
confidence level if we include the peak of low energy background events; otherwise, we face 
a 2.9 a problem. The expected energy distribution implies that the evidence for supernova 
neutrinos is at 10 a and that 1-3 background events are plausible. This conclusion does 
not change strongly when we model the time distribution of the signal, taking into account 
the presence of an initial luminous phase of neutrino emission. This suggests however that 
some of the early events (in particular, event 4) are due to supernova neutrinos and not to 
background. In summary, our comparison between the expectations and the data lead us to 
formulate the hypothesis that some of the observed low energy events are due to background 
and that some among them belong to a peculiar phase of emission, that could be further 
characterized by low energy neutrinos. Such an interpretation is mostly attractive, since it 
diminishes to a minimum the postulated number of background events and thus improves the 
agreement between the a priori (model independent) and the a posteriori (model dependent) 
expectations on the number of background events. We argue on these grounds that there is 
no significant disagreement between the average energy of the supernova neutrinos seen in 
Kamiokande-II and the conventional expectations. Alternative possibilities of interpretations 
are mentioned and briefly discussed. 
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1 Motivation and context 

On February 23, 1987, several experiments [l][2l[3j|4] contributed to begin the era of extragalactic 
neutrino astronomy. These observations had an enormous impact on astrophysics and on particle 
physics A straightforward interpretation of these data is yet difficult. This is due not only to 
the fact that we lack a firmly established theory of supernova explosion, but also to certain 
anomalous features of the data that have been understood, emphasized and analyzed in the 
course of the time. In particular, this is the case of the average energy deduced by the 12 events 
observed in Kamiokande-II (KII), which is half of that observed by 1MB, and also lower then 
what expected in theoretical models for supernova neutrino emission, recently reviewed in |14j . 
The main goal of the present paper is to propose an interpretation of this apparent discrepancy. 

1 In order to have an idea of how large was the impact it is sufficient to consider that, on March 
2007, there are 172,000 (rcsp., 3,050) entries on 'SN1987A' in Google (resp., in Google scholar). 
The SPIRES database lists 578 papers typing the same keyword; our first two references are those 
with more than 500 citations, whereas those that rank more than 100 citations include theoretical works 
on extra dimensions, CPT violation, axions, magnetic moments and exotic particles; several studies of 
the role of neutrino oscillations follow. Finally, the 20th anniversary of this observation was celebrated 
with conferences held in Moscow (Feb. 20-22, 2007), in Hawaii (Feb. 23-25, 2007) and Venice (Mar. 6-9, 
2007) j3D]. 
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The low energy feature of KII dataset is clearly reflected in the outcome of two recent analyses 
of these data, that explore certain possibilities that deviate strongly from theoretical expecta- 
tions. The first one is by Mirizzi and Raffelt [25], who describe the distribution in v e energies 
as Ep e exp[— (a + l)Ep e /(E)] and find that the best fit of KII data is provided by a ~ (i.e., 
a monotonically decreasing distribution) . The second one is the analysis of Lunardini [26] , who 
adopts a two-component distribution as suggested by three flavor oscillation scenarios and finds 
that a component with (Ep e ) ~ 5 MeV permits to fit KII data better H This type of approaches 
is useful to emphasize features of the data like the excess of low energy events in KII dataset 
but one should recall that neither (Ep e ) ~ 5 MeV nor a ~ are compatible with the current 
expectations for supernova neutrino emission, so that in a conservative analysis, these values 
should not be allowed. 

Stated otherwise, there is a trend in modern studies of SN1987A to accept the opinion that 
KII data are incompatible with the present theoretical expectations and, consequently, to start 
from this position to investigate more or less radical departures from the conventional paradigm 
for neutrino emission. One could even think that the failures to obtain a theory of supernova 
explosion could motivate such an attitude. But the fact that we do not have yet a definitive 
theory of the explosion does not mean that all theoretical possibilities are a priori equivalent. 
Furthermore, one should be careful in distinguishing between the problems of getting a theory 
of supernova explosion and the problem of knowing the distribution of the emitted neutrinos. 
In the only scenario that has been explored in some details till now, the so called "delayed 
scenario" [HI [10] , a large amount of the neutrinos-up to 90 %-is emitted in a phase that follows 
the explosion: thus, the two problems are to a certain extent independent (we will better analyze 
the connection in the following). Comparing p3] and [12] we note that the expected range of 
one crucial parameter, the average v e energies or temperatures, did not change much since 1989. 
What changed is the expectation on the temperature of the other antineutrinos and thus the 
impact of oscillations, that after [14] can be argued to modify only slightly the observable v e 
signal. To summarize, the main reasons why we are not convinced that we should abandon the 
conventional paradigm for neutrino emission are that: 1) the expectations for neutrino emission 
seem to be stable; 2) conversely, there is no convincing theoretical argument till now supporting 
an interpretation of the excess of low energy events in terms of low energy supernova neutrinos. 

Also, on general statistical ground one should be aware of the risks of using the data-12 events 
in the case of KH-to infer the characteristics of the model, rather than asking the significance 
level at which the null-hypothesis (=the theoretical expectation) is ruled out. In principle it is 
always possible to obtain a perfect fit (x 2 = 0) to the data declaring that the observed distribution 
coincides with the expected one, but only if one is ready to renounce to the previous knowledge. 
In short, we believe that a conservative discussion should address other questions: How severe is 
the deviation from the conventional theoretical expectations? How reliable is the indication of a 
large amount of low energy supernova neutrinos, inferred from KII observations? Is it possible 
to conceive other interpretations (more standard than the one proposed in [26] and [25]) of the 
observed low energy events? 

This is why we would like to explore the possibility that some of the low energy events in KII 

2 In fact, the question of which is the energy distribution of an 'average' core collapse supernova has 
important implications for future experiments: e.g., the analysis of SNf987A neutrinos of |26| was used to 
argue that future search of relic supernova neutrinos could fail because the emitted neutrinos have much 
lower energy than expected. 
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are not due to v e p — > e + n interactions of supernova neutrinos as usually assumed. The possibility 
that a few events are due to elastic scattering has been recently reconsidered [23], finding that, 
although this cannot be excluded, it helps only marginally to explain the excess of low energy 
events in KII data set. This forced the authors of [23] to admit that a few of these events are of 
a different origin, and possibly are due to background. This is a conservative position, since it 
is evident that the KII observations could be polluted by some background events: see figures 4, 
5, 10, 11, 12 of the second publication of KII collaboration on SN1987A, i.e. (5j. Furthermore, 
in Cerenkov type detector it is not possible to know whether a single event is due to signal 
or to background. In our understanding, the main point of KII collaboration was just that it 
is not possible that all 12 events are due to background, which means that an observation of 
supernova neutrinos has been made (as we show in Sect. 12. 2\ the a priori knowledge of the 
signal reinforces significantly this evidence). It should be noted that the KII collaboration did 
not quantify the probability that the 12 events are due to background on individual basis; the 
only step in this direction was to exclude occasionally one specific event, occurred below the 
threshold for solar neutrinos events (more discussion on this issue in Sect. l2~3|) . In this note, we 
analyze quantitatively the possibility that some background events occurred in the KII dataset, 
extending and complementing the analysis of [2T] . In particular, we focus our attention on the 
spatial and the energy distribution of the 12 candidate events. 

We are aware of the problems of small numbers statistics, that there are many more papers on 
SN1987A neutrinos than events, and that there is the risk of running into 'forensics' (quoting a 
witticism of John Beacom). Nevertheless, we should also recall that supernovae are rare events on 
human timescale, and that all we have at the moment is one observation of supernova neutrinos 
and some theoretical ideas to compare with. In other words, we feel that we have the duty 
to analyze all possible hints for anomalies and to extract as much information as we can from 
SN1987A neutrinos. 

The plan of this paper is the following: in Sect. [2] we discuss the tools we have to separate 
signal from background events; in Sect. [3J we formulate definite hypotheses on the background 
and on the signal events, and show that the average energy of SN1987A neutrinos agrees much 
better with the expectations when we account for the presence of a few background events; in 
Sect. H] we comment on the approach we used and the results we obtained, and finally draw our 
conclusions. 



2 Tools to separate signal from background in Kamiokande-I I dataset 



We begin recalling what are the known and the unknown aspects for the analysis of the events 
observed in KII. Several (though not all) characteristics of the background are known: the events 
are Poisson distributed with a given rate, they have relatively low energy, and they are more 
frequent in the border of the volume used in the analysis [5], that touches the planes of the 
photo-tubes. The characteristics of the signal are, instead, known theoretically: the events are 
distributed uniformly in the detector, they have relatively high energy, and their rate is not 
precisely known. Going into details, the characteristics of the background and of the signal that 
we will use are: 

1 ) The volume of the KII detector used in the analysis of SN1987A neutrinos has 



radius of R = 7.2 m, total height of H = 13.1 m, 
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for a volume of 2140 tons. This volume is much larger than the 680 tons fiducial volume defined 
in [7j, where background events are rare. We know that this wide volume is not free from back- 
ground [5] . The energy distribution of the background in our figure [2] is obtained by multiplying 
the background per unit time given in figure 2a of Lamb and Loredo [22] by the time elapsed 
from the first to the last events, 

T = 12.439 sec (2) 

that makes 2.3 background events on average and 0.272 of them above the threshold of 7.5 MeV. 
The assumption on the background has been validated by checking that the distributions of N^u 
obtained from figure 4 of [5] and the average number of background events shown in figure 10 
of the same KII paper are in good agreement with what we obtain with our distribution. More 
discussion of the assumed background is in Sect. 12.31 and App. |A] 

2) The supernova neutrino signal is supposed to be due to v e p — > ne + ('inverse beta decay') 
reaction. We use the cross section in [15] to perform the integration on the neutrino energy using 
the full matrix element and kinematical range. The antineutrino flux that we assume is the one 
described in Sects. I and II of [23], with the parameters of equation 8 there (z/„ and v T average 
energy 10 % higher than v e as in p3], energy equipartition and a = 3), including oscillations. 
The two crucial input parameters are: 

energy radiated in V e = 4 x 10 52 erg 
average energy of v e = 14 MeV 

The first value is just in the middle of theoretical expectations. Instead, the value of the average 
energy of V e is on the low side, but compatible with the present theoretical expectations [H] . 
Furthermore, it fits well the 1MB observations [2], that consisting of a sample of 8 high energy 
events are in practice background free. Two simple tests permit us to check that these values are 
reasonable: a) When the emitted antineutrino energy is compared with the (fermionic) Stefan- 
Boltzmann law Lp e = 77r 3 /240 R 2 T 4 , it corresponds to an emission from a spherical surface of 
radius R = 12 km that lasts 12.6 seconds, b) The energy radiated compares well with Lattimer- 
Yahil formulae [111 113] . In fact, assuming the total emitted energy is 6 times larger than the one 
emitted in antineutrinos (equipartition), the expected value Ej, = 1.5 x 1O 53 (M 9 /M ) 2 erg gives 
M g = 1.3M ; then the initial mass of the stellar core is M c = M g (l + O.O84M 9 /M ) = 1.4M . 
See [23] for further discussion of the selected model. With this model, we expect 11.8 signal 
events in Kamiokande-II. 

With more detailed experimental information, more precise statements could be possible; it 
would be useful to know the background rate for any single event (taking into account the specific 
position, number of hit phototubes and photoelectrons, possible correlations, etc); the error on 
the position of any event; the energy distribution of the two components of the background, 
namely the one which is located in the border (of higher energy) and the one that is distributed 
in the volume (of lower energy). With reference to the last point, we quote the Kamiokande- 
II collaboration [5]: the events with Nhu > 23 "are consistent with higher-energy products of 
radioactivity at or outside the tank wall" whereas "the events with N^u < 20 are largely due to 
214 Bi decay" (from the 222 Rn chain) and thus, presumably, more uniformly distributed. In the 
present paper we resort only to the published information on the background [3 El El [21] . 

In the following of this section, we will obtain some new hints on possible background events 
without the need of adopting a specific model for neutrino emission. However, when we will pass 
to the comparison of the expectations and the observed energy distribution, we will adopt the 
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Table 1: Value of certain observables described in the text. First line, progressive event 
number; first column, the selected observables. 



well-specified model for the signal of eq.[3]obtaining in this way other hints on possible background 
events. At this point, one could think that the conclusions that we reach on background are not 
reliable, since they depend at least partly on the assumed model for neutrino emission. Thus, 
in order to avoid confusion, we prefer to repeat that the first purpose of our analysis is to check 
whether there is a problem with the conventional expectations for the supernova neutrinos, and 
the second is to ask whether this problem can be solved by making reasonable assumptions on the 
background. So, the use of a theoretically motivated model to compare with is not only justified, 
but is in fact a necessary first step. We will stick as much as possible to the very simple model 
defined above, but will also consider reasonable (theoretically acceptable) modifications from the 
assumed model, showing that most of the results and the indications we obtain are stable. Of 
course, if a new paradigm for neutrino emission will eventually emerge, it will be important to 
update these considerations and repeat the calculations described below. 

The tools that we describe in this section are the spatial (Sect. l27Tj) . the energy distribution of 
the events (Sect. [2T2]) and the energy threshold used by Kamiokande-II collaboration (Sect. [273]) . 
The compatibility of the indications will be discussed later (Sect. [3]). 

2.1 Spatial distribution 

A lot of useful information about each of the 12 events is given in tables I and II of [5]: the number 
of hit photo-tubes Nhn, the reconstructed energy of the events (using additional information) 
and its uncertainty, the Cartesian coordinates Xi,yi,Zi (but not the error on the position) and 
finally the reconstructed direction of the momentum cos a^, cos Pi, cos 7j. In this way we can 
deduce several interesting quantities, e.g., 

6± the angle between the direction of the momentum and the normal to the closest surface of 
the volume (events with #j_ ~ are presumably seen better); 

D ou t the distance between the coordinate of the event and a formal 'exit point' from the vol- 
ume, following the direction of the momentum {Dout should be large in comparison to the 
distance between photo-tubes, about 1 meter, for an event to be reliably reconstructed); 

Di n the distance between the coordinate of the event and a formal 'entry point' in the volume, 
following again the momentum but in the opposite direction; 
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Figure 1: The volumetric distribution of the 12 KII candidate supernova neutrino events 
(dashed line) compared with a uniform distribution (continuous line). 

d m i n the minimal distance from the border of the volume used in the analysis (that can warn 
us against possible background events); 

/ the volumetric coordinate, namely the fraction of the volume contained in a concentric sub- 
detector of the same shape and with the same geometrical center of KII, that is formally 
defined as follows: 



where we set hi = Z{ — 60 cm, taking into account the offset of the origin of the coordinates. 

The last quantity is of particular interest, since when we increase the volume of the sub-detector, 
the number of supernova neutrino events should increase by the same factor (i.e., f can be con- 
sidered as the cumulative distribution coming from a uniform probability distribution function). 

In App.[B]we show that, assuming uniformity (namely ignoring systematics effects due to the 
detector response) we are underestimating the impact of a possible contamination from relatively 
high energy background events, that are known to be preferentially located in the border of the 
detector [5]. However, before proceeding we wish to stress the importance of a detailed quan- 
titative study of the detector response to a sample of uniformly distributed supernova neutrino 
events (including fluctuations, realistic light propagation, role of the energy spectrum, impre- 
cise measurement of the position, etc.) and we also warn the reader against the possibility of 
significant systematics effects, of which we are not aware. 

As can be seen in table [IJ there is no particular feature of the observable 9j_ . Similarly for the 
distance D ou t, which is usually large, except for the events 10 and 2. Anomalies emerge instead 
in the distributions of Di n , of the minimal distance from the wall d m i n and of the volumetric 
coordinate /: all these show that there are events quite close to the surface of the volume. In 
particular, d m i n is smaller than 10 centimeters for the events 2, 3, 4, 8 and 10 and / shows 
that the same events are contained in the outermost 4 % part of the 2140 ton volume, see 
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figure [TJ Performing a Smirnov-Cramer-Von Mises (SCVM) [27J [28] test for the hypothesis that 
all 12 events come from a uniform distribution we get a goodness of fit of 4.6 %. When we 
perform a test suited to reveal features in the border, this hint for a deviation from a uniform 
distribution becomes stronger. In fact, the goodness of fit is 0.3 % with the traditional version of 
the Anderson-Darling (AD) [27\ 129] test, whereas the one sided version of this test gives 0.04 %: 
see App.Ofor details. Attributing the event number 6 to background the goodness of fit becomes 
even lower, for the simple reason that this event is not located in the border. We can regard 
these results as a suggestion that there is contamination from background events located close 
(or at) the border of the volume that was used for the analysis. We interpret this as a 3 a hint, 
with reference to the numerical value of the traditional AD test. 

It is important to repeat that the assumption of a uniform distribution of neutrinos is valid 
whatever was the true (energy and/or time) distribution of SN1987A neutrinos. To the best 
of our knowledge, the analysis of the spatial distribution presented here and elaborated further 
in App. [Cjand Sect. [3] is original. In [1] we read that "the distribution of the events presented 
here is consistent with a uniform distribution", but unfortunately, the procedure adopted, the 
statistical test performed and the significance level are not indicated, so it is not possible to 
compare the results. 

2.2 Energy distribution 

In this section, we compare the observed and expected energy distributions. In the first part 
(section 12.2. ip we assess on quantitative basis the 'problem' of the energy distribution, in the 
second part (section I2.2.2j) we discuss the probabilities that the individual events to be due to 
background, that we then use to make statements on the expected number of events and on 
significance of the observations. 

2.2.1 Comparison of the expected and the measured energy distributions 

The expected event distribution is the continuous, bimodal curve in figure [2 The two distinct 
components are simply the (measured, low-energy) background and the (assumed, high-energy) 
signal due to supernova neutrinos, discussed in the beginning of Sect.[2j Clearly, the agreement of 
the theoretical distribution (the continuous curve in figure [2|) with the observations (the dashed 
curve in figure [2]) could be improved if the theoretical parameters were very different from what 
we expect, as assumed in [25] and [26]. In fact, one could be tempted to argue from figure [2] 
that it is necessary to assume a significantly lower value of the average energy of neutrinos than 
assumed in equation [3j But, is it really correct to conclude that the apparent disagreement is 
not due to a fluctuation? We performed a standard SCVM test for the null hypothesis described 
in Sect. [2] (or in figure [2]) finding that 

the assumed energy distribution (measured background + expected signal) should not 
be rejected at the 11 % significance level. 

Thus, we believe that we are justified in assuming that the theoretical distribution (that was 
motivated previously) is not in contradiction with KII observations. 

Before proceeding, however, it is useful to explain better in which sense one could claim the 
existence of a 'problem' with low energy SN1987A neutrinos. Suppose we ask whether the KII 
data could come from the assumed signal setting the background to zero] using again a SCVM 
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10 20 30 40 50 
E vis [MeV] 

Figure 2: The continuous, bimodal curve is the count rate expected in KII from the observed 
background (low energy peak) and the supernova neutrino signal, see text for a description. 
For visual comparison, we show also the spectrum obtained summing the Gaussian energy 
distributions of the 12 observed events, with and 5Ei as given in tabled (dashed curve). 



test, we find that this hypothesis can be accepted but only at 0.4 % significance level (2.9 a). 
In other words, we find that the assumed signal is in disagreement with the data if we posit 
that the background is absent. However, in view of the preceding discussion, we believe that 
this position is hardly tenable and, as we demonstrated here, it has an important impact for the 
interpretation of the data. So, we assume the concomitant presence of signal and background. 
The next step is to explore what are the quantitative implications of this position, in particular, 
for what concerns the background!! 

2.2.2 Analysis of the individual events 

In Cerenkov type detector it is not possible to state whether a single event is due to signal 
or to background. However, knowing the energy distribution of the signal and the one of the 
background, it is possible to assign to each event a probability Pi that it resulted from background 
(and a probability 1 — Pi to result from signal). 

This can be done proceeding as follows. The probability per unit second Bi (background 
rate) that an event with energy Ei results from background is known experimentally. The data 
that we use and the numerical values of the Bi are described in App. |Aj In order to obtain the 
signal rate, instead, we describe each event with a Gaussian distribution centered at the energy 

3 One may ask why we do not proceed similarly for the spatial distribution. In Sect. 12.11 we test only 
whether a subset of data is uniformly distributed (i.e., whether the subset can be attributed to supernova 
neutrino events), here we give the energy distribution of the signal plus the background and compare with 
the observations. The reason of this asymmetric procedure is simply that we do not have the a priori 
information on the spatial distribution of the background, as recalled in the beginning of Sect. [U we only 
know that the events with N^a > 23 "are consistent with higher-energy products of radioactivity at or 
outside the tank wall" and we keep in mind the message of figure 5 of [5]. 
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of the event E{ and with a width 5Ei as measured: 



Pi{E) oc exp 



2 5E? 



(5) 



Next, we calculate the convolution integral Z; of this Gaussian with the expected distribution of 
the signal. We obtain the corresponding signal rate Si = T-JT simply dividing by the time T 
from the first to the last event given in eq. O In this way we can compare the expectations for 
the background and for the signal, and we can calculate the adimensional quantity 



Pi 



B, 



Bi + Si 



(6) 



that can be thought of as the probability that an event is due to background. The approximate 
values of Pi are given below (0 means 'small' with the given precision): 



event 



12 3 

0% 1% 74% 



4 5 6 

S. 1% 93% 



7 8 9 
3% 0% 0% 



10 11 12 

51% 1% 43% 



This table indicates that the events 3, 4, 6, 10, 12 have a non-negligible chance to be due to 
background, due to their relatively low energy. 

These results do not change much when we modify the energy distribution in a manner that 
does not contradict the theoretical expectations, e.g., when we use a = 2 rather than 3, or when 
we diminish by a 1 — 2 MeV the average energy of supernova neutrinos. E.g., if the average 
neutrino energy were 12 MeV (resp., if a = 2), the convolution integral Iiq would increase by 
20 % (resp., decrease by 10 %), and thus Bio would decreases, but only by 10 % (resp., increase by 
5 %). It is pretty much evident that the probability that a low energy event is due to background 
remains quite large if the signal resembles the one shown in figure [2 

It should be noted that the number of background events above the threshold of 7.5 MeV, 
that we evaluate a posteriori using the selected model of supernova neutrino emission, is seven 
times larger than the value expected a priori. In other words, even if postulating the presence 
of a pair of background events is the best comprimise we can reach with the model, this is not 
necessarily a good compromise; this could be alternatively interpreted as an indication that the 
conventional model for neutrino emission that we adopted misses important features of the data. 
More discussion later. 

Armed with these results, we find that the probability that all events are signal is just 0.3 %, 
whereas (as obvious) it is basically impossible to have a fluctuation of 12 events that has these 
characteristics. We calculate on this basis a significance of about 10 a that at least one of these 
events are due to supernova neutrinos, that is higher than the 6.5 a significance that can be 
estimated from figure 10(b) of [5] (using only the information on the background). Though 
this result is based on the tails of the energy distributions, and should not be fully trusted 
quantitatively, it carries a reasonable and reassuring message: if we include the information on 
the expected signal, the evidence for a detection of neutrinos from SN1987A strengthens. 

The most interesting and probable cases have several background events: 



# of bkgr. events n 
probability V n 



1 2 3 4 5 
3% 26% 39% 23% 5% 
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This is easily evaluated by constructing the polynomial q[x] = IIj(l — Pj + x Pi) since, if Pi is 
the background probability, 1 — Pi can be thought of as the signal probability. The coefficient V n 
of x n counts the probability to have exactly n background events; thus we should simply expand 
the polynomial q[x] = ^2 n T ) nX n - 

Similar results are presented in table VI of [21]. But in that work the main issue is the 
study of the time distribution of neutrinos, whereas the average energy of supernova neutrinos 
is considered a free parameter whose value is decided from a global fit. Here, in order to discuss 
whether KII observations necessarily indicate an excess of low energy neutrinos, we followed a 
different approach, and adopted a (fixed, pre-selected) value of the 9 e average energy that is 
compatible with the calculated ones and with 1MB observations. Other comparisons with the 
reference results obtained by Lamb and Loredo are presented in the following. 

Let us repeat that the expectations on the background here discussed are obtained a posteriori 
adopting a pre-selected model of supernova emission. The a priori expectation on the background 
are discussed in the next Section, and a detailed comparison of the outcomes is offered in table 
[2] and discussion given there. 

In summary, we verified: (1) that the observations do not contradict seriously the theoretical 
expectations on the energy distribution of supernova neutrinos; (2) that this distribution, along 
with the measured energy distribution of the background, suggests the hypothesis that 2 to 4 
low energy events in KII dataset are due to background. 

2.3 Energy threshold 

The criterion adopted by Kamiokande-II collaboration [Tj to separate the signal from the back- 
ground is the energy threshold of 7.5 MeV (corresponding to about Nhu = 20) used in solar 
neutrino analyses. In fact, in the abstract of pQ we read: "the signal consisted of 11 electron 
events of energy 7.5 to 36 MeV"; in the text: "event number 6 has Nhu < 20 and has been 
excluded from the signal analysis". In our opinion, this criterion is perfectly fine if the aim of 
the analysis is to claim that there is an excess that cannot be explained as a fluctuation of the 
background; it is less fine if the aim is to investigate the properties of the signal. 
We propose a number of critical remarks on the exclusive use of this criterion: 

1) A priori, there is no warranty that a quantitative criterion that works well for solar neutrinos 
works also for supernova neutrino signal. In fact, solar neutrino events are directional, whereas 
u e p — > ne + events are not, the interaction rates of solar and supernova neutrinos are not the 
same, and the energy distributions are also different. However, when this criterion is adopted 
and when a time window is selected one can check the expected number of background events 
above threshold: from fig. 10 of [5], we deduce that in the time window spanned by the 12 events, 
eq. [21 this makes 0.272 events. Though this is not very large, it is not completely safe too. 

2) In general, the choice of a criterion to distinguish signal and background is not independent 
from the signal we want to reveal. For instance, the model for supernova neutrinos considered 
here suggests that the criterion of setting the threshold at 7.5 MeV is insufficient to ensure that, 
above, we have only supernova neutrinos: in fact, we find that the background is smaller than 
5 % only above 10 MeV, corresponding to about Nhit = 26 (see figure [2]). 

3) Finally, it is possible to argue against the procedure of excluding only the event number 6 
from the supernova dataset: In fact, removing the events 6, the confidence level that we have a 
uniform volumetric distribution halves when we use the AD or the SCVM test (see App. Oand 
later on). In other words, the likelihood diminishes significantly unless other signal events are 
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also assigned to background. It should be clear that this last argument is completely independent 
on the model of supernova neutrino emission. 

While we have no reservation to accept the hypothesis that the event 6 is due to background 
(see App. [A]), the arguments above suggest that it could be not the only background event in 
the time window of eq. [2[ As we discussed in the previous section, we believe that the best 
approach to extract and study the signal is to quantify the probability that an event is due to 
background; however, it is useful and instructive to consider also the simpler criterion of setting 
the energy threshold at 7.5 MeV (the comparison of these criteria will be done in Sect. 13.1.11) . 
More specifically, we will use two versions of this criterion: 

VI In the first version, we simply discard the event number 6, and consider the events that 
fall under the threshold within 1 a (namely, the events 3,4,10,12) as suspect. Thus, we will 
attempt to attribute some of them to background, in order to see if the volume distribution 
improves or not. 

V2 In the second version, we will do the same, except for assigning an additional penalty 
factor each time we attribute one of these four events to background. The penalty factor 
is simply the Poisson probability for n background events, where the average value is set 
to b = 0.272 (for instance, the penalty factor to have 2 such events is 2.8%). In this way, 
we somehow describe what we know on energy distribution of the background, though we 
are ignoring any a priori information on the energy distribution of the signal. 

These two versions of the criterion are used in table [3] of the next Section; see also the second 
line of table [2J 

3 Average energy of SN1987A neutrinos accounting for background 



As we repeatedly argued, in order to address the question on supernova neutrinos raised in the 
title of the paper we deem it necessary to discuss the role of background events in KII dataset. 
Here we collect the available information for such a discussion and draw our conclusions on 
whether there is a problem with low energy SN1987A neutrinos. As discussed in Sect. [21 we will 
assume that all the events except those that are assigned to background are due to v e p — ► ne + 
events (an extension of the 'inverse beta decay hypothesis'). Thus, the identification of the events 
due to SN1987A neutrinos becomes equivalent to the identification of the background events0 

The first task is to identify the main candidate background events, and for this aim we use 
the tools of Sect. [2j We begin listing again (a) the events that are closer to the border, namely: 
2, 3, 4, 8, 10 and (b) the events with relatively low energy: 3, 4, 6, 10, 12 with Nhit equal to 25, 
26, 16, 21, 24. Thus, we select the events 3, 4, 6, 10, 12 for further discussion: 

4 Othcr approaches arc in principle possible. For what concerns misidentification of signal against 
background, one could decide to accept only the events in the fiducial volume; if we use the definition 
of [7] we are led to keep only the five events 1, 6, 9, 11, 12. A similar possibility would be to reject all events 
under a 'fiducial threshold' that, according to our expectations, should be around 10 MeV (see Sect. 12. 3j) : 
again in this way we would remove several events including one of the previous subset. Both procedures 
should provide us with relatively safe results, but would amount to largely diminish the information on 
the characteristics of the supernova burst that we can extract from KII. In view of the small dataset, we 
believe that the most useful approach is to make an assessment on the background. 
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Events removed 


none 


6 


10 


6,10 


3,10 


3,6,10 


3,4,6,10 


(a) bkgr. above E t h 


76% 


76% 


21% 


21% 


3% 


3% 


0.3% 


(b) Vol.distr., SCVM 


4.6% 


2.4% 


12% 


7.2% 


26% 


18% 


42% 


(b) Vol.distr, AD-2side 


0.3% 


0.1% 


2.2% 


1.0% 


6.9% 


3.8% 


21% 


(b) Vol.distr, AD-lside 


.04% 


.02% 


0.4% 


0.2% 


1.5% 


0.8% 


7.8% 


(c) Energy distribution 


0.3% 


4.6% 


0.3% 


4.7% 


1.0% 


14% 


5.6% 


(d) Aver, energy [MeV] 


15 ±2.3 


15 ±2.4 


15 ±2.4 


16 ± 2.5 


16 ± 2.5 


17 ±2.6 


18 ±2.7 



Table 2: Impact of various assumptions on the background events; first row, the events 
assigned to background sample; second line, the Poisson probability of n events above 
threshold of 7. 5 MeV; third line, the SCVM significance level of a deviation from a uniform 
volumetric distribution; fourth line, the same using the two sided (traditional) AD statistic; 
fifth line, the same using the one sided (modified) AD statistic; sixth line, the probability of 
the given energy configuration, evaluated with the considerations of section \2.2\; last line, 
average visible energy of the supernova signal events. 

6 This is the lowest energy event and it has a high probability to be due to the 'diffuse' 
component of the background (e.g., internal radioactivity, radon, cosmic ray induced or 
neutrons). Being under the threshold of software analysis, Nhu = 20, it is removed in most 
investigations of supernova neutrinos (see Sect. 12.31 and App. lAj) . 

10 Besides being a low energy event, this event is known to travel a few tens of centimeters 
in the 2140 ton detector volume and to be very close to the lateral surface and to the 
upper plane of the detector-it is in the edge of the volume. We are not able to assign a 
strong quantitative significance to this information, but we believe that it puts a second 
red herring on this event. 

3, 4 These two events are quite similar; they both have low energy and are both very close to 
the border of the 2140 ton volume used in the analysis. Thus, rejecting one or the other 
from the supernova neutrino sample produces a similar effect. 

12 This is just the last event, with an energy similar to the event number 10, and thus 
(depending on the true time distribution) could fall in a region where few signal events are 
expected. It is located far from the border, though. 

3.1 A first global approach 

Now we attempt to attribute some of the events mentioned above to background, testing whether 
the spatial and the energy distributions of supernova neutrinos improve or not. The results 
are shown in table [21 where we give: (a) the Poisson probability to have a certain number 
of background events above the threshold (see Sect. 12. 3|) ; (b) the significance levels assuming 
a uniform volumetric distribution, evaluated respectively using the SCVM and the two AD 
statistics described in App. O (see also Sect. 12. ip ; (c) the probability of a given distribution of 
signal and background, that in the case where all events are declared to be signal is IT (1 — Pi) = 
0.32 %, but e.g., increases by the factor Pq/(1 — P§) = 14 removing the event number 6 (see 
Sect. I2.2[) : (d) the average energy of the signal sample, that is the aim of our discussion. 
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3.1.1 Model dependent and independent inferences 

The only really model independent test that we are aware of is the test for a uniform distribution 
of supernova neutrinos. The most powerful test between the Smirnov-Cramer-Von Mises and the 
Anderson-Darling test is the second one. In fact, this is built to reveal deviations from the 
expectations close to the boundaries, that could mean the presence of background events. We 
will use this test in the following, alone and in combination with other ones. We will also use 
the SCVM and compare the outcome of various procedures of analysis, reaching in all cases 
similar conclusions. These conclusions will be taken as a guide to formulate a hypothesis on the 
background and on the signal of KII dataset. 

Analysis A First, we present the results of the one sided (modified) Anderson-Darling test 
described in App. O With the conventional 5 % significance level, the only case among those 
of table [2] that should not be rejected is the one where all four candidate background events are 
attributed to background (see fifth row of table [2]). Another interesting conclusion that we obtain 
from the one sided AD test is that at least 3 events close to the border should be attributed 
to background. The conclusions till here are independent on the energy distribution. When we 
use the a priori expectations on the energy distribution of the background we conclude that the 
lowest energy events, among those selected by the AD test, have a higher chance to be due to 
background. This is the closest we can go to a model independent argument. 

Analyses B1,B2,B3 Next, we quantify the improvements (likelihood ratios) in the description 
of the supernova neutrino signal under the various hypotheses on the background described in 
Sect. El We will compare four different procedures. In all of them, we perform a traditional 
(two sided) AD test to check the uniformity of the assumed supernova events in the volume. 
Additionally: 

Bl) we follow Kamiokande-II and assign a priori the event number 6 to background (see Sect. 12.31 
first version of the threshold criterion); 

B2) additionally to B2, we impose a Poisson penalty for background events above the threshold 

of 7.5 MeV (see Sect. 12. 3| second version of the threshold criterion); 

B3) we test the likelihood of the given energy distribution as in Sect. 12.21 

With any of the four procedures we evaluate the factors of improvements of the given hypothesis, 
choosing as a comparison the case when all 12 events (or 11 events in the second case) are 
considered to be due to supernova neutrinos. E.g., the procedure B3 suggests that the hypothesis 
that only the event 6 is due to background is more probable than the hypothesis that all events 
are signal, by a factor of improvement equal to 

P„ ;.(all events but 6) x P ere .(all events but 6) ^ ^ 
P vo i (all events) x P en .(all events) 

(approximate values of P are given in table [2]) . The results of these calculations are given in 
table O 

The procedure Bl makes very little use of the background and of the signal; the procedure 
B2 uses the a priori knowledge on the background, but neglects any information on the signal; 
the procedure B3 uses all we know, though its reliability depends also on the model for supernova 
neutrinos we assume. It is difficult to make firm conclusions from this table, but we note that 
even in the more pessimistic procedure (B2) the presence of several background events cannot 
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AD-2side + 


6 


10 


6,10 


3,10 


3,6,10 


3,4,6,10 


remarks 


Sect. [2,31 VI 


= 1 




8.2 




30 


160 


model indcp. 


Sect.[UV2 


EE 1 




2.2 




1.1 


0.5 


model indcp. 


Energy distrib. 


5.7 


7.0 


48 


64 


500 


1100 


model dependent 



Table 3: Factors of improvement in the description of the Kamiokande-II supernova neu- 
trino events for the selected hypotheses on the background. The subset of events that are 
assigned to background is in first line. In addition to the Anderson-Darling test for a 
uniform volumetric distribution, we use the information listed in the first column. When 
we impose the Kamiokande-II threshold (second and third rows) the event 6 is discarded a 
priori from the supernova dataset. Last column, comments on the selected procedure. 



be firmly excluded; this is because an a priori unlikely fluctuation of the number of background 
events above threshold is compensated by the fact that some of these events lie close to the 
border of the detector. On the other hand, we see that with the procedure of analysis B3, 
the only one that makes a full use of the knowledge on the background and on the signal, the 
factors of improvements are about a thousand times when all 4 candidate background events are 
assumed to be background, and half of that value when we assign events 3, 6, 10 (or 4, 6, 10) 
to background. In short, table [3] suggest that a compromise between the extreme assumptions — 
namely, of considering completely unknown or reasonably known the energy distribution of the 
signal-would admit the presence of several background events. It is important to note that the 
results of the more aggressive procedure of analysis (B3) do not change much when we modify 
the energy distribution in a manner that does not contradict the theoretical expectations, see 
Sect. 12.21 In other words, one could be tempted to argue that the model dependence of our 
stronger analysis is not large, if the astrophysical parameters lie in the range suggested by the 
theory and by the 1MB observations: see again Sect. [2] and 



Analyses C1,C2 Finally, one may ask what happens when we apply the SCVM test for the 
spatial distribution of the events, instead than the AD test that we used so far, or when we 
change the procedure of analysis. We show that the conclusion is very similar to the previous 
ones, though it is true that the weight of the energy distribution becomes comparably more 
important with the SCVM test. 

Let us begin with a description of the procedure CI. According to our hypothesis, any of the 
12 events can be due either to background or to inverse beta decay signal. Thus we consider the 
2 12 = 4096 individual possibilities, and for any of them we test the likelihood of the given spatial 
distribution, using this time the SCVM test and the likelihood of the given energy configuration 
as we did in Sect. 12.21 (we checked that with a SCVM test we obtain similar outcomes). Next, we 
consider the product of the two likelihoods and normalize the sum of the likelihoods of the 4096 
mutually exclusive cases to unity. We find that the occurrence of several background events is 
likely: 



# of bkgr. events n 
probability 



1 2 3 4 5 6 

0.1% 2.2% 14.6% 37.3% 33.8% 11.4% 0.5% 



The individual cases with probability above 5 % are the following ones: 
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Volume 



Energy 



Angle 




0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 



Figure 3: Cumulative distributions for the case when the events number 3, 4, 6 are assigned 
to background, and significance of the SCVM test. In the abscissa the value Fj of the 
cumulative distribution (see App. lC\) . 



background events 


3,6,10 


3,4,6,10 


3,4,6,10,12 


3,6,10,12 


3,4,6 


3,6 


probability V n 


19.1% 


18.3% 


9.9% 


8.1 % 


7.7 % 


6.8 % 



We note that the event number 12 is occasionally assigned to background because of its low 
energy. What about the event number 6? In all those cases that have probability larger than 3 
%, this event is due to background. We also remark that the 2048 cases when this event is due 
to background have a total probability of 90%. 

In the procedure of analysis C2 we test also the SCVM likelihood of the angular distribution 
following [23] . The events that are most commonly assigned to background are the usual ones, 
namely those figured out in the beginning of Sect. [3l The individual cases with probability above 
5 % (the most probable individual cases) are those when the background is given by the events 
number 3,4,6 (15.4 %), 3,6 (14.0 %), 3,6,10 (11.9 %); 3,4,6,10 (9.1 %), where in brackets we give 
again their a posteriori probabilities, that by definition are normalized to give 1 when we sum 
over the 4096 possible cases. The cumulative distributions in the volume, in the energy and in 
angle for the most probable single case (the first one listed just above) are shown in figure El 
along with the likelihood of the three experimental distributions. Remarks on the persistent 
discrepancy of the angular distribution will be offered in Sect. [H The overall probabilities to 
have 2, 3 or 4 background events becomes 28%, 41% and 20% and the probability that all events 
are due to signal is 0.5%. 



3.1.2 A tentative explanation of low energy events 

In view of the arguments above discussed, we formulate our hypothesis on the signal and on the 
background: 

9 (or less) among the 12 events observed in Kamiokande-II were due to SN1987A 
neutrinos, whereas the remaining 3 (or more) low energy events were due to back- 
ground. 

The motivations and the merits of this position are discussed in Sect. 13.1.11 

Quite importantly, this hypothesis allows us to argue that there is no excess of low energy 
supernova neutrinos. Consider as a quantitative indicator the average visible energy of super- 
nova neutrino events, that being weighted with the detection efficiency and the cross section is 
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significantly higher of the average antineutrino energy-see, e.g., [23]. With the signal described 
in the beginning of Sect [2J the average visible energy of supernova neutrino events is expected 
to be 



where n is the number of events due to supernova neutrinos. If we assume that all events are due 
to signal (n = 12), we find that an average energy that is 2.6 a lower than the expectations. This 
discrepancy is very similar to the 2.9 a 'problem' of low energy neutrinos that was formulated in 
the beginning of Sect. 12. 2t one can say that the selected observables catches the crucial feature 
of the 'problem'. Next we deduce the value of the average visible energy of supernova neutrino 
events from the data, in the hypothesis on the signal formulated just above. In the case when 
we assign the events 6, 10, 3 (resp., 6, 10, 3, 4) to the background the average visible energy 
of supernova neutrino events is 1.4 (resp., 0.9) standard deviations below the expected value 
(approximate values are in the last line of table [2]). In other words, they are in reasonable 
agreement. Thus, if 3 or 4 low energy events were caused by background the question raised in 
the title of the paper should receive a negative answer. 

3.2 A second global approach 

The previous interpretation leaves something to be desired. In fact, we expect only 0.272 back- 
ground events above the threshold of 7.5 MeV on average in the time between the first and the 
last event, so there was only a 3 % probability a priori that 2 or more such background events 
occurred (see third row of table [2]). In other words, even if the invoked background fluctuation 
cannot be just dismissed, it is significantly larger than a priori expected. Furthermore, the event 
number 10 is in a region of the detector much noisier than on average (compare table [H Sect. [3] 
and fig. 5 of [5]) but the events 3 and 4 seem to be located in a relatively safer region. In 
other words, the hypothesis discussed previously is the best one within the considered theoretical 
context, but it is not particularly good: this can be seen quite clearly from figure [21 that shows 
that the low energy events we are discussing are just in between the signal and the background 
distributions. So, we are motivated to consider whether we can evade the conclusion we reached 
and/or how we can reach a more satisfactory conclusion. We will show that in order to address 
at least partially these needs it is sufficient to consider a more accurate model for the supernova 
neutrinos, in particular, for what regards the time evolution of the neutrino luminosity. This 
possibility is quite appealing, since it is largely within the conventional expectations. 

3.2.1 A more refined hypothesis on neutrino luminosity 

A phase of intense neutrino luminosity is expected during the first hundred of milliseconds from 
the collapse. This phase has a marked non-thermal character [8] and it is thought of as to 
play a key role for the explosion [U [10] , [18] . In the following we term this hypothetical phase 
as "accretion", and the subsequent thermal phase as "cooling". The latter phase is the one 
considered till now. For the cooling phase, we will continue to assume equipartition, the same 



E vis = 20.6 MeV. 



(8) 



This can be compared with the value given from the data, using the formula 




(9) 
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Figure 4: Cumulative time distribution for the 16 events seen by KII in a time window of 30 
seconds, compared with the expected distributions for an emission with constant luminosity 
(long dashed) and with a 'standard' accretion phase (continuous line). 



amount of emitted energy, the same average energy as in eq. El and an emission time of the order 
of the time from the first to the last event (T in eq. [2]). Accretion instead lasts just a fraction 
of second and is much more luminous. We further assume that during accretion equipartition is 
violated by an excess of electronic neutrinos and antineutrinos: to be definite, we suppose that, 
on average, there are 3 v e and 2 v e each Vu->v T ,Vu,, or u T . We describe the v e luminosity in the 
simplest manner we can conceive: 

{— for t < r a ("accretion" or non-thermal phase) 

I »=t , (10) 

— x e T c for t > r a ("cooling" or thermal phase) 

where we select reasonable values for the parameters, namely £ a = 8 x 10 51 erg and r a = 0.5 s for 
accretion, and £ c = 4x 10 52 erg and r c = T for cooling. With these numerical values, we can derive 
several other interesting quantities: 1) the total emitted energy £j, = 4.5£ a + 6£ c = 2.76 x 10 53 
erg; 2) the fraction emitted during accretion A.5£ a /£b = 13 %; 3) the electronic lepton number 
emitted during accretion, that corresponds to about half solar mass of iron if the average energy 
of u e is of 12 MeV; 4) the energy of 20 times 10 51 erg (=20 foe=20 bethe) emitted in v e and P e , 
that could be sufficient for the explosion. None of these numbers should be taken too literally, 
but they are comfortably within the present theoretical expectations. 

Let us come now to the interpretation of the events. Since the average energy is supposed to 
stay unchanged and the total emitted energy is almost the same, the total number of supernova 
events and their energy distribution do not change much: only their time distribution changes. 
Before explaining why this position is relevant to the interpretation of the low energy KII events, 
we begin noticing that the expected rate of event accumulation, 



N(< t) = 0.187 t + 



2.0 — for t < T a 



2.0 + 11.8 1 -exp 



T a ~ t 



(11) 

for t > T a 
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compares reasonably well with the KII data, since the (SCVM) GOF is 7.7%. In the previous 
equation we show the contributions of the background-first term-and of the signal separately; the 
number of events during cooling is the same as before, 11.8. For comparison, the hypothesis that 
the signal accumulation rate is 11.8/T for t < T-that is, supernova neutrinos are emitted with a 
constant luminosity in an interval of T seconds-has a GOF of 1.9% only. See figure HI where we 
show the sequence of events observed by Kamiokande-II in a time window of 30 seconds. The last 
events are well explained by the expected background curve; instead, the remarkable feature of 
data close to t = is only partially explained by our assumption that a standard accretion phase 
of neutrino emission also occurred. Indeed, the expected number of events during accretion is 
only 2.0 (a bit less than the value 2.4 = 11.8 x £ a /£ c due to oscillations). If it was possible 
to double the number of events during accretion, the GOF would increase to 30%. The other 
parameters of emission play a less crucial role for the subsequent discussion^ Similarly, if we 
would arbitrarily omit the event number 3 and the event number 6 from the dataset, the GOF 
would increase to 25%. For other recent analyses of time distribution that also find an evidence 
of an initial luminous phase see [20], EH [24]. Now, we pass to discuss the impact of the updated 
hypothesis on the interpretation of the events. 



3.2.2 Accretion and the meaning of the early events 



The theoretical context defined in section [3.2.11 suggests strongly the possibility that some of the 
early events seen in KII occurred during a phase much more luminous than average: £ a /T a ~ 
5 x £ c /t c . Identifying, arbitrarily but not unreasonably, the arrival of the first neutrino to the 
Earth with the time of first detected event, this possibility concerns the events number 1,2,3,4. 
Thus, the analysis of section 12.21 should be updated, taking into account that the signal rate 
in the denominator depends explicitly from the time t{ at which the individual event occurred. 
More precisely, the factor Si in equation [6] should be multiplied by £(U) that depends on the 
instantaneous luminosity: 

' 4.15 



for t < T a 



cxp 



for t > T a 



(12) 



Of course, the value of S a /r a (the luminosity during accretion) determines the numerical factor 
4.15, and the value of the luminosity in the beginning of the cooling phase £ c /t c gives = 1 
just after t = r a . A straightforward calculation gives the new values for the probability Pj that 
an individual event is due to background: 



event 
Pi 



1 



3 
41% 



4 5 6 
3% 1% 94% 



7 



10 11 12 

67% 2% 67% 



(compare with Sect. l2.2T2l) . The most striking consequences of the new hypothesis is the increased 
probability that the event number 4 is due to supernova signal, whereas the late events have a 
higher chance to be caused by background. Due to these features, the expectations on the number 
of background events, based on the energy distribution, are not radically modified: 



# of bkgr. events n 
probability V n 



1 2 3 4 5 
7% 28% 41% 20% 2% 



5 It is possible to somewhat improve the fit to the data using slighly different model parameters (e.g., 
a value of r c some times smaller than we assumed) but this would not change radically our conclusions. 
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The probability that all events are signal is 0.36%, and we can exclude at 10 a that all events 
are due to background. However, there is an important difference with what we found based 
on the results of section 12.21 The events 3,4 (10) that are located close to the border of the 
detector have a smaller (larger) chance to be due to background, and this has an impact on 
global analyses. Using as a reference analysis C2 described above, we find the following results, 
that can be compared with those in section 13.1.11 

• The probability that no background events occurred becomes 0.9%. This is small is abso- 
lute terms, but larger than estimated above using only the energy distribution. 

• The probability that a certain number of background events occurred is: 



# of bkgr. events n 
probability V n 



1 

11% 



4 5 
14% 2% 



that is rather similar to the values calculated just above, based on the expected energy 
distribution only. 

The individual cases that occurred with probability above 5% are: 



background events 


6,10 


3,6,10 


6,10,12 


6 


3,6 


3,6,10,12 


probability V n 


18.9% 


12.0% 


9.5% 


7.9 % 


7.0 % 


5.7 % 


volumetric SCVM 


7.2% 


18.0% 


3.3% 


2.4 % 


6.5 % 


9.9 % 


angular SCVM 


3.3% 


1.2% 


1.8% 


8.4 % 


4.0 % 


5.2 % 



where we give also the GOF of the volume and of the angular distributions in the last two 
lines. Note that these 6 cases alone cover 61% of the 4096 possibilities. 



The event number 6 is due to background in 



of the cases. 



In summary, the existence of a accretion phase brings a new twist in the analysis of the back- 
grounds, and leads to somewhat different conclusions especially regarding the interpretation of 
the event number 4. However, the new analysis still leads us to expect that some low energy 
events among the events number 3,6,10,12 could be attributed to background. 



3.2.3 Again on the interpretation of low energy events 

In view of the discussion of the previous section, the interpretation proposed in section [3 . 1 . 21 does 
not need a real revision, except, possibly, replacing the event 4 with the event 12 in the role of 
candidate background event. 

However, the wider theoretical context also allows us to take a new point of view, or in other 
words, it is possible to formulate another way out from the 'problem' of low energy events. In 
fact, one could argue that the events occurred during the non-thermal phase (accretion) should 
not be compared with the ones occurred later. In particular, the average energy of the first four 
KII events is 12.6 ±2.4 MeV, that is pretty much lower than the one of the next seven events (the 
event number 6 being attributed to background), namely 17.1 ± 3.3 MeV. Of course this is due 
to the presence of a few low energy events among the first events (the events 3 and 4). But while 
in section 13.1.21 we regarded them as as possible background events, here we are lead to suggest 
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that they (or at least some of them) are due to a peculiar phase of neutrino emission, that in 
principle could be characterized by a lower average energy than during thermal phase. In short, 
the fact that KII measured several low energy events can be at least partly attributed to the fact 
that some of these events occurred during accretion. In order to have an acceptable volumetric 
distribution, it is sufficient to assign only one additional event between 3 and 10 (along with 
event 6) to background. 

4 Discussion 

We have shown in Sect. [3] that the agreement between the observed and expected average energy 
of supernova neutrinos improves significantly assuming that 2-4 of the low energy events in KII 
are due to background, and that this interpretation is supported by the volumetric distribution 
of the events. A few of these low energy events could also belong to a peculiar phase of neutrino 
emission ("accretion"), contributing to solve the 'problem'. In other words, we believe that the 
null hypothesis (i.e., the minimal model) is not significantly challenged by the excess of low 
energy events observed by Kamiokande-II (or: we have enough freedom to evade the conclusion 
that there is a problem with the conventional expectations). Here, we would like to offer various 
remarks on this conclusion. 

New tools used. A new crucial information that permitted us to reach new conclusions is the 
analysis of the spatial distribution of the events (Sect. 12 . X [) . The validity of the hypothesis on 
the background that we formulated in Sects. I37L21 and [3.2.31 rests also, to a certain extent, on the 
assumption we made on the signal of supernova neutrinos described in the beginning of Sects. [2] 
and 13.2. IT In our view, these should be considered a conservative and reasonable assumptions on 
the signal. They are consistent with the 1MB observations that are basically background free. 
Furthermore, the conclusions we reach remain valid under small modifications of the model for 
supernova neutrinos, and do not contradict (but rather strengthen) what we obtain following 
Kamiokande-II and setting the energy cut at 7.5 MeV (Sect. f2~3T> . 

Peculiarity of the approach. The simple but crucial observable that we have used to draw our 
conclusions is the average value of the visible energy, see in particular Sects. 13.1.21 and 13.2.31 A 
detailed comparison of the energy distributions in figure [2] would suggest two anomalies, namely: 
a) a defect of high energy events around 20 — 40 MeV and b) an excess of events around 10 MeV. 
While we argued that the low energy anomaly could be due to background events (or to the fact 
that several events occurred during accretion) we are unable to explain the high energy anomaly 
if not resorting to a fluctuation of the data. However, 'global' observables are safer against the 
effects of the fluctuations; thus, in view of the low number of collected events we emphasized the 
comparison of average visible energies rather than the comparison of energy distributions. In 
this respect, our formulation of the 'problem' of low energy events, or equivalently our approach 
to KII dataset, are more conservative and distinct from the ones of other recent analyses such 
as |25] and [26], and closer to the one of [23] (see figure 1 there). A quantitative statement is in 
the first paragraph of Sect. 12.2.11 

Leftover possibilities and/or alternative (but not exclusive) keys of interpretation. It is not 
impossible that other reasonable effects contribute to explain the apparent excess of low energy 
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events (and/or defect of high energy events) in KII dataset, for instance: 

(i) The measured energies should not be thought of as Gaussians (as suggested by the reported 
data) . If we use errors that scale as y/~E [19] rather than being constant as assumed in the rest of 
the paper (equation [5]) we see a modest improvement of the agreement in the region dominated 
by the background (i.e., low energies). 

(ii) There is one event due to elastic scattering that degrades the visible energy as allowed by 
the selected theoretical model of emission [23]. This option is interesting since when assigning 
events to background, the angular distribution does not improve much (see e.g., figure [3]). 

(Hi) The neutrino flux in the thermal phase deviates slightly from expectations, say, in the 
direction suggested by [25] but possibly not that radically (in our evaluations, the change is not 
large-see Sect. I3.1.1|) . 

(iv) Some observed event is really due to supernova neutrinos, but it suffers of a very poor 
energy reconstruction (though, we do not have the necessary information to elaborate on this 
possibility) . 

(v) Some observed events are due to supernova neutrinos, not necessarily u e , interacting with 
the walls of the detector (a similar possibility was originally noticed for non-standard collapses 
|16| , but was shown to yield a relevant amount of events in existing detectors also for a standard 
collapse [T7]). 

(vi) The time structure of the signal is significantly different from what we know, e.g., the phase 
of accretion is much more luminous (see Sect. 13.2. 1 [) . 

While most possibilities are largely within conventional expectations, some of them go beyond 
the minimal model that we formulated. This is true for the last two possibilities listed just above, 
that are however tightly connected with the most unknown aspects regarding the explosion, and 
(if proven to be viable) could offer us new insights on what happened during SN1987A neutrino 
burst. 

Possible future developments. With more information, more detailed analyses of the KII data 
could be possible. In fact, one should analyze the energy and the spatial distribution at the same 
time, assuming that the observations are due to several phenomena; a signal of high energy due 
to neutrino interactions with free protons, with electrons and with oxygen nuclei, uniformly dis- 
tributed in the detector (and possibly, also a hypothetical component of the signal concentrated 
on the wall of the detector); a low energy background component, similarly distributed, with (in 
principle) known intensity; a more energetic background component, concentrated in the wall 
of the detector, also known (again, in principle). More discussion of these needs is in App. lAl 
and App. [HI Another interesting thing would be the calculation of the best fit for the supernova 
neutrinos parameters, keeping into account the presence of background in KII, with peculiar 
(and in principle known) spatial and energy distributions. For lack of information-as discussed 
in the beginning of Sect. [2j compare also with the discussion of equation 3.22 in |21|-we have 
been forced to analyze the energy and the spatial distributions separately, simply distinguishing 
in each analysis between a signal and a background component. Perhaps, more detailed analyses 
like these will be carried out in the future; in the meantime and view of what we learned here, we 
are lead to expect as a plausible outcome that the event 6 should be attributed to background 
because of its low energy, and some events among 3, 4, 10 (resp., 10, 12) should be attributed to 
background mostly because of their spatial distribution (resp., of their relatively late occurrence). 

We wish to close recalling the warning of E N Alexeev: "it is possible that some of the low 
energy events carry an important message, still to be understood". We believe that these are 
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1.6 5 


1.9 3 


2.9 2 


1.2 2 


2.1 3 


3.7 2 


4.5 5 


8.2 5 


1.5 5 


1.5a 


1.9 3 


I.62 


3.8 2 


2.9 2 


2.8 2 


3.8 2 


B* 


1.0 5 


5.4 4 


3.1 a 


8.5 3 


5.3 4 


7.1 a 


5.0 6 


1.0 5 


1.0 5 


1.8a 


4.O4 


1.4a 


7.3 2 


5.2 2 


I.82 


7.3 2 



Table 4: 1 st line, the progressive event number; next, the estimated background rate in 
Hz/MeV according to the procedure of the 'convolution integrals' followed by LL (2 nd line) 
and the procedure we describe in the text (3 rd line). A subscript indicates the exponent: 
e.g., 32 means 3 • 10~ 2 . 

wise words and we subscribe them; as we saw, the message could be the existence of a luminous 
phase of neutrino emission (Sect. 13.20 . However, we also believe that it is important to check 
how close (or how far) we can go to describe the data using only conventional hypotheses on 
supernova neutrinos, and this is what we tried to do in this paper. Certainly, the handful of 
events collected from SN1987A was enough for the first observation of supernova neutrinos (the 
quantitative statement on the significance of KII observation is in Sect. 12.23]) . but as we showed 
with our analysis, the possibility that some of these events are not due to v e p — > n e + should 
suggest caution whenever we attempt to infer the characteristics of SN1987A neutrino emission 
from the observations. 
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A The assumed background 

We interpret the curve B(E) presented in fig. 2 of |21] as the measured background rate for the 
estimated event energy. We tested this interpretation by constructing a histogram from fig. 4 
in [5] and using a linear conversion factor between N^u and E, fixed by the correspondence 
between known values: N hit = 20 44> 7.5 MeV and/or N hit = 26 44> 10 MeV. Within limited 
statistics, the histogram coincides with the curve given by Lamb and Loredo (LL). Note that 
our interpretation agrees with the one of LL [21], who state "KII and Baksan provided us 
with measurements of B(E) n and later describe their fig. 2 as "background rate measurements". 
Now, since the measured background already includes the effects of the fluctuations in the energy 
measurements ('smearing'), we estimate the background rate for the individual event using the 
central value of the measured energy: Bi = B(Ei), see table SJ A different prescription was 
followed in [21], namely, B\ L was evaluated by taking the convolution integral of the curve B(E) 
with the Gaussian energy distribution of the individual events; the values from table III of [21] 
are given for comparison. In our understanding, such a prescription leads to double counting 
the effect of 'smearing', and thus (1) to underestimate the B{ for the events with energy close to 
the background peak, and (2) to overestimate those in the tails of the B(E) distribution]! The 

6 We believe that the procedure of taking the convolution integrals should be applied to the true energy 
distribution of the background, namely, the distribution that does not include the fluctuations in the energy 
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new values of Bi have been obtained as follows: below 10 MeV, we obtain accurate data from 
the plot of LL, whereas, above, this is not of much use; above 9 MeV, however, we can use the 
background curves denoted as STEP1 and STEP2 shown in fig. 3 of [6] for the purpose of solar 
neutrino studies, with the conversion between energy and N^a given there (STEP1 agrees with 
the curve of LL, where comparable). Below 7 MeV, the new procedure leads to higher Bi (ev.s 
6, 13-16); from 7.5 to about 9 MeV, there is no significant difference with the procedure of LL 
(ev.s 3, 10, 12); between 9 ad 17 MeV, the background curves in [6] suggest a lower Bi (ev.s 
2,4,5,11); above-in the most uncertain but safest from the background point of view, and also 
less crucial region for parameter estimation-we assume a lower background rate than assumed by 
LL (ev.s 1, 7-9). More precise values could be provided by the experimental collaboration using 
B(E, x)-not the average value B(E)— thus describing the role of the position of the individual 
events (similarly for detection efficiencies). It should be noted that the changes in the values of 
the Bi that we propose are not crucial for the analysis; their main effect is just to increase the a 
priori probability that the event number 6 is due to background, as in the first and most popular 
interpretation [I]. 

B The assumption of uniform volumetric distribution 

We analyze critically the hypothesis that the events from the signal are uniformly distributed. 
In the relevant publications [IJ [5j [7J [6] we were not able to find warnings against important 
systematics of this type, but we tentatively identified three such physical effects that could affect 
the observed signal events: 

Geometry: The events occurring close to the border and propagating toward the closer wall could 
be missed, since the Cerenkov light could pass through the phototubes without being recorded. 
In our understanding, this is the main reason why the detection efficiency of Kamiokande-II 
never reaches 100%, even for the higher energy events, but is at most 92%. In order to test 
this hypothesis, we considered the external part of the detector that has a distance from the 
phototubes smaller than their typical distance of 1 meter. Since this part has a volume of about 
800 m 3 , we reproduce the 92% efficiency if ~ 20% of the events are lost in this way, that is a 
reasonable figure being comfortably smaller than 40%. Again, this effect leads us to expect that 
the loss of signal events is more significant for the events produced close to the border, especially 
for / > 0.6, that corresponds to the external part of the detector considered just above. We 
attempt to quantify the size of this geometrical effect by the function G[f] = 1— (/— 0.4) 2 #[/— 0.4], 
where we used a quadratic function to describe a smooth behavior {9[x\ is the Heaviside function). 
This means that ~ 8% of the total number of events are systematically lost close to the border, 
and that even at / = 1, the loss is smaller than 40%. 

Light attenuation: The light from the low energy events can be lost because of light attenuation, 
where A ~ 50 m (in [5] , we read that A "exceeds 50 m at any time" ) . Let us evaluate the average 
attenuation coefficient for Kamiokande-II detector, as a function of the volumetric coordinate /. 
Considering the simplified case of propagation in straight rays, we find the differential distribution 
in the variable / by averaging over the other coordinates that specify the position of the event 

measurement. The convolution integral would then amount to include the effect of the fluctuations on 
individual basis, and could lead to a more accurate background evaluation. 
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and the direction of emission: 



d 4 = /'* f% /' dl ,^["/ 1/3 ./ 1/3 .^]+^[/ 1/3 ."/ 1/3 .^] (13) 

df JO JO Jo 37T 

where A[r,h,c,cp] = (e~ Din / x + e~ Dout ^ x )/2 is the average attenuation for a given direction of 
observation, Di n and D ou t are the two distances from the walls of the detector (see also table [1]), 
and we measure the radial coordinate r in units of R and the height h in units of H/2. A 
straightforward calculation gives the approximation for dA/df ~ 0.858 + 0.068 /, that means 
that the events produced in the center suffer an attenuation effect that is ~ 10% larger. The 
systematic effect due to light attenuation depends on the energy. If we assume that the signal 
is uniformly distributed in energy, we are lead to expect that the loss of low energy events from 
the center is more significant. 

Spectrum: However, the previous point does not take into account the real spectrum of the 
signal. In order to demonstrate that this is also relevant, it is enough to consider a simple one 
dimensional situation, when the light propagates just in one direction or the opposite one. Let 
us compare what happens to the events produced in the middle (/ ~ 0) and in the border 
of the detector (/ ~ 1). In the first case, the light attention effect is e~ L / 2X ~ 0.87, where 
L = 14 m is the linear size of the detector; in the second case, this depends on whether the event 
propagates toward the closer border (no attenuation) or toward the opposite border, when we 
have e~ L l x « 0.76. In the first case, we lose all the events between N m i n and 1.15A^ m j n , where 
Nmin is the minimal number of photons that gives a trigger; in the second one, we lose half of 
the events between N m i n and 1.3iV m j n (where e L l x « 1.3). If the signal is uniformly distributed 
in energy, we see that at 0{L/\) the effect is identical in the two cases; thus, the effect described 
by eq. [13] is purely three dimensional. But if the distribution of the signal increases with the 
energy, as we expect, the loss of signal events is more significant for the events produced close 
to the border. In fact, the events that fall below the detection threshold in the first case (events 
produced close to / ~ 0) are those that belong to the least populated part of the spectrum. E.g., 
if we assume that the energy spectrum of the signal rises linearly from zero, the total number of 
remaining events in the region N m i n — 1.2>N m i n is 50% larger in the first case. This last effect is 
presumably the less relevant one. 

In summary, various effects can change the distribution of the detected signal events: the 
second effect increases a bit the number of expected supernova neutrino events close to the border, 
the other two diminish it. It is not possible to estimate precisely the systematics deviations from 
uniformity without a detailed simulation of the detector response that takes into account the 
energy distribution of the signal, the role of energy fluctuations, etc. However, if we account 
for the two systematics effects that work in opposite directions simply by multipling G[f] by 
dA/df[f], we see that the geometrical effect is the most important one. This is the reason why 
we claim that the GOF with a uniform volume distribution underestimates a possible problem 
with events close to the border: in fact, when we apply the SCVM test to the 12 events seen 
by Kamiokande-II we get a GOF of 3.4% rather than the value that we obtain with uniform 
distribution, 4.6% (see table [2]). Similarly, when we discard event 6, the value of the GOF 
becomes 1.8% rather than 2.4% (again in table [2]). As we see, the changes due to a more 
accurate modeling of the detector response are not expected to be very large, and are in any case 
smaller than the changes due to the use of alternative statistical tests (see next Appendix). 
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C The Anderson-Darling test 



The Anderson-Darling (AD) test is used to test if a sample of data comes from a specific 
distribution. It is a modification of the Kolmogorov-Smirnov test and gives more weight to the 
boundary region than the Smirnov-Cramer-Von Mises test. 

Consider a set of data x±,X2, arranged in increasing order. Given the expected cumu- 
lative distribution function F[x], we calculate the values in the points xf 



Fj = F[ Xj ] 



(14) 



Next, consider the empirical cumulative distribution function -F n [a;], that counts the number of 
events below x; thus, F n [x] = j/n for Fj < u < Fj + i, where we set Fq = and F n+ \ = 1. The 
quantitative indicator of the likelihood of the given data proposed by Anderson and Darling is: 



W 2 =n (F n [x] 
Jo 



u) 2 ip[u] du 



(15) 



F\x] 



where ip is a positive function on the interval (0, 1). The case when ip = 1 gives the well-known 
Smirnov-Cramer-Von Mises test, that somehow resembles the x 2 statistics: 



1 



scvm 



3=1 



3 ~ 1/2 



n 



(16) 



Let us pass to the alternative possibilities. The traditional (two sided) version of the Anderson- 
Darling test statistic is obtained setting ip[u] = l/(u(l — u)), that gives: 



traditional 



-n 
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£ 

3=1 



2j-l 



n 



\og[F 5 (1 - F n+ x-j)] 



(17) 



This provides a particular sensitivity close to the points u = and u = 1. A trivial modification 
consists in setting ip[u] = 1/(1 — u): in this way, the test becomes particularly sensitive close to 
u = 1. The modified (or one sided) Anderson-Darling test statistic can be written as: 



modified 



2j-l 



3=1 



n 



lOgfl - Fj 



n+l-j\ 



(18) 



Since l/(u(l — u)) = 1/u + 1/(1 — u), equation [T71 can be recovered by summing a third test 
statistic when we replace ip[u] = 1/u, that would be suited to test for features close to u = 



1 



n+l-j 
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(the practical formula for this third test statistic is obtained replacing Fj 
equation \TE\\ . 

A null-hypothesis that we want to test is whether a pre-selected subset of the KII data is 
uniformly distributed in the detector, as expected if they are due to supernova neutrino events. 
Thus, in our case the number of data is n < 12, the coordinate x is what is called / in Sect. 12. 1\ 
and the cumulative distribution function is simply F\x\ = x. The AD tests of Eqs. [T7] and [TH] 
are of particular interest since the points close to / = 1 are in the border of the volume used in 
the analysis, and, there, it is not possible to exclude a priori an effect of contamination from the 
background. 
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Given the test statistics, we calculate the value W 2 for the subset of KII data in which we are 
interested. How often this value results from a uniform distribution? This can be obtained with 
a straightforward procedure: 1) we extract randomly n data from a uniform distribution for iV 
times (N = 10 8 in our case); 2) we calculate the value of W 2 for each data set; 3) finally, we test 
the frequency of the condition W 2 > W 2 , namely, the number of times N* that this condition is 
satisfied over the number of extractions N. 

E.g., we infer from the test of equation [T7] that when all 12 events (resp., all except event 
6) are supposed to be due to supernova neutrinos, we get W 2 = 4.95 (resp., W 2 = 5.80). A 
value this high can occur as a result of a statistical fluctuations, but only in 0.32 % (resp., 0.13 
%) of the cases; this means that we can reject the hypothesis that the 12 (resp., 11) events are 
uniformly distributed at the 2.9 a (resp., 3.2 a) significance level. Similarly, we infer from the 
test of equation [T8l that the case when all events (resp., all except event 6) are assumed to be 
signal can be rejected at the 3.5 a (resp., 3.8 a) significance level. 

In principle, it is possible to object that the situations when different statistical tests give 
very different outcomes should be regarded as dubious and that the statistical criterion would be 
better selected a priori rather than a posteriori; these are the reasons why in this paper we do 
not use exclusively the (more powerful) AD test, but also the (more conventional) SCVM test. 
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